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Abstract 

We have investigated the lattice thermal transport across the asymmetric tilt 
grain boundary between armchair and zigzag graphene by nonequilibrium molec- 
ular dynamics (NEMD). We have observed significant temperature drop and 
ultra-low temperature-dependent thermal boundary resistance. More impor- 
tantly, we find an unexpected thermal rectification phenomenon. The thermal 
conductivity and Kapitza conductance is direction-dependent. The effect of 
thermal rectification could be amplified by increasing the difference of temper- 
ature imposed on two sides. Our results propose a promising kind of thermal 
rectifier and phonon diodes based on polycrystalline graphene without delicate 
manipulation of the atomic structure. 
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Graphene, a two dimensional honeycomb lattice, is well recognized as one of the 
most promising materials in future electronics and nanotechnology[Tl [2]. Nu- 
merous previous studies have shown that the monolayer graphene samples with 
perfect lattice structure exhibit outstanding electronic, thermal, mechanical and 
optical properties 131 SI [3 IS]- While the size of graphene flakes from mechani- 
cal exfoliated is too small to be utilized for mass production of functional de- 
vices, chemical vapor deposition (CVD) enables the synthesis of large-area, high- 
quality and low-cost graphene sheetslJllHlE]- Due to the imperfections of metal 
foils in the growth process, the grains inevitably nucleate during the growth 
which results to the formation of grain boundaries in graphene sheets [TUl [TT] . 
Understanding the effect of the grain boundaries on the mechanical, electronic, 
and thermal properties is very helpful and necessary for future applications of 
the larger-scale polycrystalline graphene. 

Recently, several studies have shown that the grain boundaries have cru- 
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cial influence on the mechanical and electronic properties. For the mechanical 
property, the grain boundaries with high density of topological defects show 
anomalous characteristics of mechanical strength[T2]. In the case of charge 
transportation, the charge mobility depends on the grain size which is deter- 
mined by the grain boundaries pi. Furthermore, the transmission of the charge 
carrier across the grain boundary depends on the symmetry of the atomic struc- 
ture. The charge transport of the symmetric grain boundary reveals high trans- 
parency while that of the asymmetric grain boundary emerges almost perfect 
reflection[13]. Experiments have shown that the line defects like grain bound- 
ary can be one of the most powerful ways to manipulate electronic properties at 
nanoscale [141 [T51 [TB] . For the thermal transport, although there have been many 
studies on the point defects including vacancies [T7], Stone- Wales defects [TB] and 
isotope defectsfW, the effect of line defects like grain boundary on thermal con- 
ductivity is still far from well-understood. In a recent work j20], thermal trans- 
port across the tilt grain boundaries with simplest symmetric atomic structure 
has been considered. Sudden temperature drop at the grain boundary and su- 
perior boundary conductance has been observed in that system[20j. There is 
still a long way for us to further explore and understand the impact of one- 
dimensional global defects on the thermal transport of graphene. Similar to the 
charge transport, it is reasonable to expect that the asymmetric grain boundary 
would exhibit different property of thermal transport from that of the symmetric 
grain boundary. 

In this paper, we adopt the nonequilibrium molecular dynamics (NEMD) 
method[2Tl HH US] to simulate the heat transport across the asymmetric grain 
boundary between armchair graphene and zigzag graphene. We show that the 
Kapitza conductance and the thermal conductivity have strong dependence on 
the temperature. Furthermore, surprising thermal rectiflcation phenomenon 
has been observed in the graphene with such asymmetric grain boundary. Our 
results provide a practical choice for designing thermal rectifier in polycrystalline 
graphene without delicate atomic engineering. 

Figure 1 shows the atomic structure of the asymmetric grain boundary 
connecting the armchair and the zigzag regions. This structure has been de- 
noted as a boundary structure between the armchair graphene and the zigzag 
graphene jUj while the electronic transport across that has been studied [1^125]. 
This structure is composed of a periodic array of 5-pentagon and 7-heptagon 
topological defects. According to the mismatch of the period of the zigzag 
part and the armchair part, this asymmetric grain boundary shows peculiar 
7-5-7 structure while the standard pentagon-heptagon dislocation cores would 
fiip their orientations from 5-7 to 7-5[24 . In spite of the lattice mismatch, the 
formation energy of this asymmetric grain boundary is much smaller than the 
typical energy of the bare edges in graphene [26J. 

In all simulations, velocity Verlet algorithm is used to integrate the equation 
of motions [27]. Periodic boundary condition is imposed on the direction parallel 
to the grain boundary. Along the direction perpendicular to the grain boundary, 
fixed boundary condition is applied on the green atoms in the Figure 1. The 
brown atoms in Figure 1 are coupled to Nose — Hoover heat baths f28], whose 
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temperatures are set to Thigh and Tiow , respectively, to generate the temperature 
gradient. The length of heat bath used in our simulation is around 1 
nm. We have carefully tested the effect of the heat bath length on our 
results. The quantity of the calculated thermal rectification is only 
slightly affected by the length of heat bath. A recent parameter-modified 
Tersoff potential is used to describe the atomic interaction [21], which has been 
verified to give the right velocities of the acoustic-phonon which is crucial in 
thermal transport. The new parametrized potential considerably improves the 
computational values of the lattice thermal conductivity and makes them in 
good agreement with the experimental results [50] . After the system reaches 
to the steady state, the heat flux could be calculated by the energy flowing from 
the heat bath to the systempi]. The thermal conductivity could be obtained 
according to the Fourier's law 

(1) 



ATiwh 

In the equation, J represents the heat flux, d is the length of sample, ATi 
is the temperature difference between Thigh and Tiowi w and h is the width 
and height of the monolayer graphene, respectively. The temperature usually 
shows significant temperature jump at the grain boundary region, the boundary 
conductance is defined as: 

«=Ak <^> 

J is the heat fiux and AT2 is the temperature jump at the boundary|20|. 

Before starting the simulation, the atomic structure has been fully optimized 
in a microcanonical ensemble by conjugate gradient method for ensuring the 
forces on each atom being less than 10~^ eV/A. Applying the Nose — Hoover 
thermostat on two ends about 1 ns, the system reaches steady equilibrium. The 
heat fiux and the temperature gradient is obtained by the average of 4 ns. The 
time step of the simulation is set to 1 fs. The simulation time has been tested 
to be enough for convergence. The width of the period parrallel to the grain 
boundary is about 3 nm and the height of the sample is set to be the distance 
between graphite layers about 3.35 A [20]. 

At first, we calculate the thermal conductivity of pure graphene to validate 
our simulation results. While the length of the sample is about 10 nm and 
25 nm under the simulation temperature around 300 K, the thermal conduc- 
tivity is about 327 W/(mK) and 455 W/(mK) respectively, which is in good 
agreement with the previous NEMD simulation results [32] [33] [31] . Under the 
same simulation condition, the thermal conductivity of the graphene sheet with 
the asymmetric grain boundary reduced to only around 240 W/(mK) which 
becomes considerably lower than the pure graphene. The apparent reduction 
of thermal conductivity is attributed to the strong phonon scattering at grain 
boundary. 

Superior boundary thermal conductance has been observed in the simula- 
tion. At room temperature, the boundary thermal conductance of the asym- 
metric grain boundary is about 1.5 x lO^*' W/m^K which is much larger than 
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the reported value of grain boundaries in nanocrystalline diamond thin films 1 35] 
and silicon-silhcon (001) S29 grain boundaries [5^ 157]. Our results is slightly 
lower than the previous reported boundary thermal conductance of the sym- 
metric grain boundarv|20j. The misorientation angle of this asymmetric grain 
boundary is 30 degree which is larger than that of the reported symmetric grain 
boundaries. The large misorientation angle would induce the increase of the 
defects densities and the probability of the phonon-boundary scattering, which 
would finally reduce the thermal boundary conductance. 

The temperature dependence of the boundary thermal conductance and the 
thermal conductivity has been investigated systematically, which are shown in 
Figure 2(a) and 2(b), respectively. When the temperature increases up to 550 K, 
the thermal boundary conductance increases from 1.3x 10^° W/m^K to 2. Ox 10^° 
W/m^K. The similar trend has also been observed in the Si-Ge interfaces [55] 
and the Kr-Ar interfaces |39]. While at the high temperature around 600K, the 
boundary thermal conductance suddenly drops. This could be attributed to the 
transformation of atomic structure of the grain boundary. The C-C bonds in 
the array of 5-pentagon and 7-heptagon topological defects broken has been ob- 
served which reveals the structural instability of the asymmetric grain boundary 
at high temperature. The thermal transport depending on the atomic structure 
of the grain boundary was also found in a previous study |2D]. The thermal 
conductivity follows the same relation as the boundary thermal conductance 
respecting to the temperature, which has been reported in previous theoretical 
study on defect free graphene |17Llin] . while the sudden drop of the thermal con- 
ductivity near 600 K could be also attributed to the atomic structural changing 
of the grain boundary. 

The monotonic increase of thermal boundary conductance with the rising of 
temperature before the structure changing could be interpreted by considering 
the following two factors. Firstly, the thermal fluctuation increases with the rise 
in temperature. Consequently, the probability of inelastic phonon scattering 
increases at the grain boundariesjSS]. The optical phonon with high frequency 
might break down into numerous acoustic phonons with low frequency. These 
acoustic phonons have higher probability of transmission comparing to high 
frequency optical phonon due to the limited width of the asymmetric grain 
boundary. Besides that, more acoustic phonon would be stimulated at high 
temperature which could also contribute to the increase of boundary thermal 
conductance. 

It is interesting that we find the thermal conductivity and Kapitza con- 
ductance is asymmetric, i.e., the heat fiux from armchair to zigzag is larger 
than that from zigzag to armchair, which suggests that such asymmetric grain 
boundary structure could be used to design a new kind of thermal rectifier. 

The thermal rectification is found to be strongly temperature dependent and 
it could reach to as high as 74%. The thermal rectification versus the tempera- 
ture difference |Ar| and the average temperature T are plotted in Figure 3(a) 
and 3(b). In Figure 3(a), the temperature difference between Thigh and Tiow 
has been set to 120 K and the simulation temperature varies from 200 K to 700 
K. With the decrease of the simulation temperature, the thermal rectification 



5 



increases from 1.05 to 1.73. In Figure 3(b), the average simulation temperature 
iom+ hzgh "i^^g been set to 400 K and the temperature difference between Thigh 
and Tiow varies from 40 K to 140 K. With the increasing of the temperature 
difference, the thermal rectification increase from 1.05 to 1.33. In both Figure 
3(a) and 3(b), the thermal rectification increases with the rise in the temper- 
ature difference. The nanostructure like the carbon nanotube intramolecular 
junctionjni], the asymmetric graphene ribbon jU], the carbon nanocone|42| and 
the hybrid graphene-graphane nanoribbon[43] has been proposed to be thermal 
rectifier in the past few years. The effect of thermal rectification in the asym- 
metric grain boundary revealed in this work is comparable to those in previous 
proposed systems. 

Furthermore, the thermal rectification is actually determined by 
the change of the phonon spectral overlap under opposite bias tem- 
perature[46]. In fact, the obvious difference between the phonon 
power spectrum of the armchair and zigzag graphene was verified in 
the previous theoretical studies [47, 48]. Besides that, the temperature- 
dependent phonon power spectrum, a typical characteristic of the 
nonlinear system, would change the phonon spectral overlap while 
the temperature bias is reversed, which results in the emerging of 
thermal rectification. On the one hand, if the armchair part is under the 
high temperature, the power spectrum between these two parts is well matched 
and the thermal transport could be efficiently exchanged between two parts. On 
the other hand, if the armchair part is under the low temperature, the power 
spectrum of the armchair and the zigzag parts are shifted and the overlap is ex- 
pected to be diminished. By the decrease of overlapping phonon spectrum, the 
efficiency of the thermal exchange between the armchair part and the zigzag part 
declines dramatically. While enlarging the temperature difference in two ends, 
the mismatch between the phonon spectrum of two parts would also increase 
therefore the thermal rectification is magnified. Under low temperature, the 
thermal conductivity of the graphene highly depends on the temperature |40| . 
which indicates that the phonon spectrum shifted more considerably at low tem- 
perature than at high temperature by the same amount change of temperature. 
For this reason, with the same temperature difference, the overlap of phonon 
spectrum decreases more at low temperature than that at high temperature, 
therefore the effect of thermal rectification is amplified at low temperature. 
According to the phonon spectrum match/mismatch theory mentioned above, 
the results of our simulation could be well understood. On the other hand, 
at the high temperature region more acoustic phonons would gen- 
erate in the system. If more acoustic phonons are generated, the 
grain boundary with limited width would impose less impact on the 
thermal transport. In this case, the thermal rectification would also 
decrease at the high temperature region. 

In summary, we have found that the Kapitza conductance across the grain 
boundary and the thermal conductivity are strongly temperature-dependent. 
Surprisingly, the asymmetric heat flux across the asymmetric zigzag-armchair 
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grain boundary is observed. It provides a potential candidate for future design 
of carbon-based thermal rectifier which has attracted a lot of attention recently 
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Figure 1: Structure of the asymmetric tilt grain boundary used in simulation. The green atoms 
are set to be applied fix boundaries. The brown atoms are set to be coupled to heat bath. The 
yellow atoms are the asymmetric boundary structure composed of 5-pentagon and 7-heptagon 
topological defects. The size of the sample is defined by the length between two heat bath. 
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Figure 2: The temperature dependence of (a) the Kapitza conductance and (b) the thermal 
conductivity of armchair-zigzag asymmetric grain boundary. Both the Kapitza conductance 
and the thermal conductivity increase with the temperature increasing. The sudden drop of 
the Kapitza conductance and the thermal conductivity is attributed to the changing of the 
boundary structure at high temperature. In all the simulations, the size of the sample is 10 
nm. 
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Figure 3: The thermal rectification versus (a) different temperature and (b) different temper- 
ature difference. In Figure 3(a) the temperature difference is set to be 120 K. In Figure 3(b) 
the simulation temperature is set to be 400 K. The effect of thermal rectification has been 
evaluated by the ratio of the asymmetric heat flux along two reverse directions. means the 
heat flux from the armchair part to the zigzag part wrhile J_ means that the heat flux from 
the zigzag part to the armchair part. In all the simulations, the size of the sample is 10 nm. 



